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ABSTRACT  For  split  Hopkinson  pressure  bar  (SHPB)  tests  on  soft  materials,  the  goals  of  homogeneous  defor¬ 
mation  and  uniform  uniaxial  stress  in  the  specimen  present  experimental  challenges,  particularly  at  higher  strain 
rates.  It  has  been  known  for  some  time  that  attainment  of  these  conditions  is  facilitated  by  reducing  the  thickness 
of  the  specimen  or  by  appropriately  shaping  the  loading  pulse.  Typically,  both  methods  must  be  employed.  Pulse 
shapes  are  often  tailored  to  deliver  a  smooth  and  sufficiently  slow  rise  to  a  constant  axial  strain  rate,  as  this  promotes 
equality  of  the  mean  axial  stress  on  the  two  faces  of  the  specimen,  a  condition  referred  to  as  dynamic  equilibrium. 
However,  a  constant  axial  strain  rate  does  not  eliminate  radial  acceleration,  which  may  result  in  large  radial  and 
hoop  stresses  and  large  radial  variations  in  the  radial,  hoop  and  axial  stresses.  An  approximate  analysis  (assuming 
homogeneous  deformation  and  incompressibility)  indicates  that  these  radial  inertia  effects  would  be  eliminated  if  the 
radial  strain  rate  were  constant.  Motivated  by  this  result,  we  consider  loading  pulses  that  deliver  a  constant  radial 
strain  rate  after  an  initial  ramp-up.  The  corresponding  axial  strain  rate  is  no  longer  constant  on  any  time  interval, 
but  for  sufficiently  thin  specimens  the  resulting  departure  from  dynamic  equilibrium  may  be  small  enough  to  be 
tolerable.  This  is  explored  here  by  comparing  the  analytical  predictions  for  the  conventional  and  “optimal”  loading 
pulse  shapes  with  corresponding  numerical  simulations  of  SHPB  tests  on  a  soft,  nearly  incompressible  material. 

1.  Introduction 

The  split  Hopkinson  pressure  bar  (SHPB),  also  known  as  the  Kolsky  bar,  is  widely  used  to  characterize  the  strain  rate 
sensitivity  of  inelastic  materials  in  a  state  of  compressive  uniaxial  stress  for  moderate  to  large  strains.  The  review 
articles  by  Gray  [1],  Gama  et  al.  [2],  and  Ramesh  [3]  and  the  recent  book  by  Chen  and  Song  [4]  are  recommended 
for  background  on  this  technique.  The  standard  specimen  shape  for  SHPB  tests  is  solid  disk,  which  (particularly  for 
soft  specimens)  is  relatively  thin  compared  to  the  diameter  in  order  to  reduce  axial  inertia  effects.  The  analysis  and 
simulations  in  this  paper  are  confined  to  this  specimen  geometry.  The  standard  analysis  of  SHPB  test  data  relies  on 
the  assumption  of  uniform  uniaxial  stress  and  uniform  biaxial  strain  throughout  the  specimen.  These  conditions,  as 
well  as  a  nearly  constant  nominal  strain  rate,  can  often  be  achieved  after  an  initial  ring-up.  However,  a  number  of 
technical  challenges  arise  when  the  specimen  is  extremely  soft;  cf.  [1]  [4]  and  also  Gray  and  Blumenthal  [5],  Chen  et 
al.  [6],  Song  and  Chen  [7],  Moy  et  al.  [8],  Song  et  al.  [9],  and  Sanborn  [10]. 

Inertial  effects  in  soft  specimens  may  result  in  highly  non-uniform  conditions  and  a  stress  state  that  is  far  from 
uniaxial,  particularly  if  the  strain  rate  is  sufficiently  high.  To  minimize  the  effects  of  axial  inertia,  loading  pulse 
shapes  are  typically  tailored  to  deliver  a  smooth  and  sufficiently  slow  rise  to  a  constant  nominal  (or  engineering) 
axial  strain  rate.  This  promotes  axial  uniformity  of  the  stress  and  strain  components  and,  in  particular,  equality  of 
the  mean  axial  stress  on  the  two  faces  of  the  specimen,  a  condition  referred  to  as  dynamic  equilibrium.  The  constant 
axial  strain  rate  condition  also  simplifies  calibration  of  constitutive  models  for  which  the  stress  depends  explicitly 
on  the  strain  rate,  e.g.,  plasticity  models  with  rate  dependent  yield  stress.  However,  a  constant  nominal  axial  strain 
rate  does  not  eliminate  the  radial  acceleration  of  the  specimen,  which  may  result  in  large  radial  and  hoop  stresses 
and  large  radial  variations  in  the  radial,  hoop  and  axial  stresses.  Since  only  the  mean  axial  stress  is  measured  in  an 
SHPB  test,  the  full  stress  state  cannot  be  inferred  in  such  cases;  in  particular,  the  deviatoric  stress  (which  typically 
exhibits  the  largest  rate  effects)  cannot  be  determined.  Consequently,  the  measured  axial  stress  is  not  useful  for 
constitutive  model  calibration  unless  the  data  can  be  “corrected”  for  radial  inertia  effects. 

An  approximate  analysis  (assuming  an  incompressible  specimen  that  deforms  homogeneously)  indicates  that  the 
radial  inertia  effects  would  be  eliminated  if  the  nominal  radial  strain  rate  were  constant.  Motivated  by  this  result,  we 


consider  loading  pulses  that  deliver  a  constant  nominal  radial  strain  rate  after  an  initial  ramp-up.  The  corresponding 
axial  strain  rate  is  no  longer  constant  on  any  time  interval;  in  fact,  it  is  necessarily  decreasing  on  the  time  interval 
over  which  the  radial  strain  rate  is  constant.  But  for  sufficiently  thin  specimens  the  resulting  departure  from  dynamic 
equilibrium  may  be  small  enough  to  be  tolerable.  This  is  explored  here  by  comparing  the  analytical  predictions  for 
the  conventional  and  “optimal”  loading  pulse  shapes  with  corresponding  numerical  simulations  of  SHPB  tests.  An 
inertial  correction  for  the  deviatoric  stress  is  also  compared  with  the  deviatoric  stress  in  the  simulation. 

The  analysis  on  which  these  the  inertial  corrections  and  the  optimal  loading  pulse  shapes  are  based  applies  only 
to  nearly  incompressible  specimens,  that  is,  materials  for  which  the  bulk  modulus  is  several  orders  of  magnitude 
larger  than  the  shear  modulus.  This  is  not  necessarily  a  drawback,  since  radial  inertia  effects  are  expected  to  be 
more  pronounced  for  nearly  incompressible  specimens.  Examples  of  soft,  nearly  incompressible  materials  include 
many  rubbers,  biological  materials  with  high  fluid  content  (in  particular,  brain  tissue),  and  tissue  surrogates  such 
as  gelatins.  The  analysis  assumes  that  the  specimen  is  isotropic,  although  no  particular  form  for  the  constitutive 
relation  is  required.  On  the  other  hand,  an  isotropic,  nonlinear  elastic  constitutive  relation  was  used  for  the  specimen 
in  the  numerical  simulations,  namely,  a  compressible  version  of  the  Mooney-Rivlin  model  with  the  bulk  modulus  of 
water  and  a  shear  modulus  typical  of  a  gelatin,  so  that  the  ratio  of  bulk  to  shear  modulus  exceeds  104. 

2.  Homogeneous,  Axisymmetric  Deformations 

Our  analysis  is  confined  to  the  deformation  and  stress  in  the  specimen,  which  is  assumed  to  be  a  (relatively  thin) 
solid  disc.  We  use  a  cylindrical  coordinate  system  aligned  with  the  axis  of  symmetry  of  the  specimen  and  the  pressure 
bars.  Coordinates  of  the  specimen  in  the  undeformed  reference  configuration  are  denoted  by  (R,  O,  Z);  coordinates 
of  the  corresponding  material  point  in  the  deformed  state  are  denoted  by  (r,  0,  z).  For  a  general  axisymmetric 
deformation, 

r  =  f(R,Z,t),  0  =  0,  z  =  z{R,Z,t),  (1) 

where  t  denotes  time.  The  radial,  hoop,  axial  and  shear  components  of  the  Cauchy  or  true  stress  tensor  <x  (taken 
positive  in  compression)  are  denoted  by  oyr,  agg,  ozz  and  <rrz.  The  other  two  stress  components,  arg  and  <jzg, 
are  zero  for  axisymmetric  deformations  since  the  specimen  is  assumed  to  be  isotropic,  but  in  general  (e.g.,  if  the 
specimen  bulges)  arz  need  not  be  zero.  The  Cauchy  stress  measures  force  per  unit  deformed  area.  The  nominal 
(engineering,  1st  Piola-Kirchhoff)  stress  tensor  £  measures  force  per  unit  original  (undeformed)  area.  The  axial  and 
radial  components  of  this  tensor  (which  are  the  only  ones  referred  to  in  the  sequel)  are  denoted  by  T,zz  and  Er^  and 
are  likewise  taken  positive  in  compression. 

We  make  the  approximation  that  the  specimen  deformation  is  also  homogeneous,  in  which  case  (1)  reduces  to 

r  =  Xr(t)R,  0  =  0,  z  =  Xz(t) Z  +  z0(t) .  (2) 

The  initial  radius  of  the  specimen  is  denoted  by  Rq,  and  the  deformed  radius  at  time  t  is  denoted  by  ro(t): 

r0(t)  =  Xr(t)R0  and  r/r0(t)  =  R/Rq  .  (3) 

For  a  material  point  located  initially  at  radius  R ,  the  deformed  radius  r  at  time  t  is  given  by  (2)i,  so  the  radial 
velocity  and  acceleration  of  the  specimen  at  this  radial  location  are  given  by 


i~'  =  :''(t)K=yfTr- 


where  a  superposed  dot  denotes  a  time  derivative.  At  any  instant  t,  the  radial  velocity  and  acceleration  of  the 
specimen  increase  linearly  with  the  radius.  They  are  zero  on  the  axis  and  attain  their  largest  absolute  value  at  the 
(stress-free)  lateral  surface,  where  r0(t)  =  Xr(t)Ro  and  r 0(t)  =  Xr(t)Ro .  From  the  relations  (2),  we  have 


dr  r 

Xr  =  M  =  R=:Xe 


l  \  _  dz 

cinci  Xz  —  • 

oZ 


It  follows  that  Ar,  Xg  and  Az  are  the  stretches  (local  ratios  of  deformed  to  undeformed  length)  in  the  coordinate 
directions.  Since  we  are  concerned  with  possibly  large  specimen  deformations,  we  must  be  careful  to  distinguish 
between  various  finite  deformation  measures  of  strain.  The  two  most  commonly  used  measures  in  the  Hopkinson 


bar  literature  are  the  nominal  (or  engineering)  strain,  which  is  the  change  in  length  per  unit  original  length,  and  the 
logarithmic  (or  true)  strain.  The  nominal  and  logarithmic  measures  of  radial,  hoop  and  axial  strain  are  defined  by 

er  =  l  Ar ,  eg  =  1  Xg ,  ez  —  1  \z ,  (6) 


and 


£r  =  —  In  Ar  ,  £0  =  -lnA0,  £2  =  -lnA2, 

respectively.  Since  the  radial  and  hoop  stretches  are  equal,  so  are  the  corresponding  strains: 

eg  =  er  and  eg  =  er  ; 


(7) 

(8) 


so  we  will  focus  on  the  radial  and  axial  stretches  and  strains.  The  nominal  and  logarithmic  strain  measures  approach 
each  other  in  the  small  strain  limit,  but  differ  substantially  for  large  strains.  We  will  work  primarily  with  the  nominal 
strains  in  this  paper,  as  this  simplifies  most  of  the  results.  Note  that  stretches  are  1  in  the  undeformed  state,  greater 
than  1  in  extension,  and  between  0  and  1  in  compression.  The  strains  have  been  defined  so  that  they  are  positive  if 
the  particular  coordinate  direction  is  in  compression  and  negative  if  it  is  in  extension.1  In  an  SHPB  test  the  specimen 
is  compressed  in  the  axial  direction,  so  that  0  <  A2  <  1  and  0  <  e2  <  1.  Since  the  lateral  surface  is  unconfined,  the 
specimen  is  free  to  expand  radially  and  we  expect  that  Ar  >  1  and  hence  er  <  0.  From  (6)  and  (7),  we  see  that  the 
axial  strain  rates  and  stretch  rates  are  related  by 


Analogous  relations  hold  for  the  radial  strain  and  stretch  rates. 


(9) 


Now  consider  the  stress  state.  Let  p  denote  the  pressure  and  s  the  deviatoric  part  of  the  Cauchy  stress  tensor.  Then 


(Jrr  —  Srr  T  p  ,  &6Q  —  $66  T  Pi  & zz  —  &zz  T  P  •> 

Srr  T  §66  T  Szz  —  0  ,  p  —  3  ( &  zz  T  (Jrr  T  (JQO^  • 


(10) 

(ii) 


Since  the  radial  and  hoop  stretches  coincide  and  the  specimen  is  isotropic,  the  radial  and  hoop  stresses  are  equal  and 
the  shear  stress  is  zero: 

arr  =  agg ,  oy2  =  0 .  (12) 

Thus  a  homogeneous  axisymmetric  deformation  results  in  a  biaxial  stress  and  strain  state;  the  principal  axes  of  stress 
and  strain  are  the  axis  of  symmetry  and  any  axes  orthogonal  to  it.  In  this  case,  the  relations  (11)  simplify  to 

Srr  —  Sgg  2  $zz  5  V  —  3  (Tzz  T  2 (Jrr )  •  (18) 

On  substituting  the  relation  (13)2  for  p  into  the  relation  (10)3  for  cr22,  we  obtain  the  equivalent  relations 

<?zz  =  I  Szz  +  (Trr  ,  SZZ  =  §  (<TZZ  ~  CTrr)  ■  (14) 

For  axisymmetric  deformations  of  an  isotropic  specimen,  the  balance  of  radial  momentum  is  given  by2 


dar 


(Jrr  —  crgg 


dor 


=  -pr , 


dr  r  dz 

where  p  is  the  density  in  the  deformed  state.  In  view  of  (12),  for  homogeneous  deformations  this  reduces  to 


(15) 


d(Jrr 

dr 


—p  r . 


(16) 


Since  the  lateral  surface  of  the  specimen  is  stress  free,  arr  =  0  when  r  =  ro(t).  Then  on  integrating  (16)  from  an 
arbitrary  radius  r  to  ro(t)  and  using  (4)  and  (3),  we  obtain 


2  (T (t ) 


2  ( 7 j'f  (^)  1 


(17) 


1The  standard  sign  convention  in  the  continuum  mechanics  literature  takes  stress  and  strain  components  as  positive  in  tension.  The 
opposite  sign  convention  is  typically  used  in  the  compression  Hopkinson  bar  literature,  and  we  have  followed  that  convention  here. 

2 The  minus  sign  on  the  right  side  of  (15)  is  a  consequence  of  the  convention  that  stresses  are  positive  in  compression. 


where 


(18) 


arr(t)  =  Xr(t)Xr(t) . 

Working  directly  with  (17),  we  find  that  arr(t)  is  the  mean  value  of  arr  over  the  deformed  volume  of  the  specimen 
at  time  t.  Since  the  deformation  is  homogeneous,  this  is  also  the  mean  value  of  arr  over  the  undeformed  reference 
configuration.  And  since  arr  is  axially  uniform,  arr  is  also  the  mean  value  of  arr  over  any  deformed  (z  =  constant)  or 
undeformed  (Z  =  constant)  cross-sectional  area.  Since  agg  =  arr ,  we  conclude  that  the  radial  and  hoop  stresses  vary 
quadratically  with  the  radius  but  are  axially  uniform.  They  attain  a  peak  absolute  value  of  2arr  ( t )  at  the  centerline 
(r  =  R  =  0)  and  reduce  to  zero  at  the  lateral  surface.  They  have  the  same  sign  as  the  radial  stretch  acceleration 

Ar,  but  without  additional  assumptions  neither  the  sign  nor  the  magnitude  of  Xr  can  be  inferred.  In  any  case,  it 
is  clear  that  the  presence  of  nonzero  radial  and  hoop  stresses  is  a  radial  inertial  effect,  that  is,  a  consequence  of  the 
fact  that  the  pr  term  in  (16)  is  not  necessarily  negligible.  It  follows  that  in  the  quasi-static  limit,  that  is,  in  the 
limit  as  r  (or  Xr)  approaches  zero,  arr(t)  =  0.  Hence,  as  expected,  in  a  quasi-static  test  we  have  a  uniaxial  stress 
state  (oyr  =  agg  =  0),  and  by  (13)2  and  (14)2  the  pressure  and  deviatoric  stress  are  completely  determined  by  the 
axial  stress:  p  =  \ozz  and  szz  =  | azz.  On  the  other  hand,  it  is  clear  from  (13)2  and  (14)2  that  these  quasi-static 
estimates  for  p  and  szz  will  be  in  error  if  arr  is  sufficiently  large  relative  to  azz. 

3.  The  Incompressibility  Approximation 


The  Jacobian  of  the  deformation,  denoted  by  J,  is  the  determinant  of  the  deformation  gradient  F  and  represents  the 
local  ratio  of  deformed  to  undeformed  volume.  For  a  general  deformation,  J  is  the  product  of  the  principal  stretches; 
for  a  homogeneous  axisymmetric  deformation  this  yields  J  =  XrXgXz  =  A,2AZ.  Since  the  focus  of  this  paper  is  on 
nearly  incompressible  specimens,  we  will  make  the  approximation  that  the  deformation  is  volume-preserving.  Then 
J  =  1,  which  is  equivalent  to  any  of  the  following  relations: 
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Thus  the  radial  stretches  or  strains  are  determined  by  axial  stretches  or  strains,  and  vice  versa.  In  this  case  the  true 
and  nominal  stress  components  are  related  by3 


o\zz  =  XzTizz  =  (1  —  ez)T,zz  ,  Orr  =  Ar£r_R  .  (20) 

Now  recall  the  relations  (9)  between  the  axial  components  of  strain  rate  and  stretch  rate  (with  analogous  relations 
for  the  radial  components).  By  (19)  we  obtain  the  following  additional  relations  between  the  radial  and  axial  rates: 
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From  the  relations  on  the  left  in  (21),  we  find  that  the  radial  strain  and  stretch  accelerations  are  given  in  terms  of 
the  axial  stretch  or  strain  rates  by 
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We  also  have 
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although  this  relation  is  less  useful  than  (22) .  For  later  use  in  the  discussion  of  optimal  pulse  shapes,  we  note  that 
the  nominal  axial  strain  acceleration  is  given  in  terms  of  the  radial  stretch  rates  by 


V  o  c  7'r) 
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(24) 


Now  consider  the  radial  stress.  Since  the  deformation  is  volume-preserving,  p  =  p0,  the  density  in  the  undeformed 
state.  On  substituting  this  and  the  relations  (19)i  and  (22)  for  Ar  and  Ar  into  (18),  we  obtain  the  following 


3These  relations  follow  from  the  general  relation  a  =  J  1'SFT ,  or  from  the  relations  between  deformed  and  undeformed  area. 


Figure  1:  Axial  strain  and  radial  stretch  rates  (a)  and  accelerations  (b)  for  a  smooth  ramp-up  to  a  constant  nominal 
axial  strain  rate  4  of  2500/s  after  135  /is.  Also  shown  in  (b)  is  the  normalized  mean  radial  stress  (dashed  line). 


expressions  for  the  mean  radial  (and  hoop)  stress  in  an  isotropic  specimen  undergoing  a  homogeneous,  volume¬ 
preserving,  axisymmetric  deformation: 
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The  radial  (and  hoop)  stress  distribution  is  then  given  by  (17).  For  a  given  radial  or  axial  strain  history,  we  see 
that  the  mean  radial  stress  at  any  instant  is  proportional  to  p0Rq .  Thus  radial  inertial  effects  can  be  reduced  by 
decreasing  the  radius  of  the  specimen.4  Note  that  the  only  material  property  appearing  in  (25)  is  the  density  p0 ;  in 
particular,  this  estimate  for  the  radial  stress  is  independent  of  the  constitutive  relation  for  the  specimen.  The  relation 
(17),  with  arr(t)  given  by  the  expression  on  the  right  in  (25),  is  equivalent  to  relations  in  Dharan  and  Hauser  [11], 
Warren  and  Forrestal  [12],  and  Scheidler  and  Kraft  [13]. 


Now  consider  a  conventional  smooth  loading  pulse  for  an  SHPB  test  with  rise  time  t\.  We  take  time  t  =  0  to 
be  the  instant  at  which  the  loading  pulse  arrives  at  the  specimen-incident  bar  interface,  so  that  e2(0)  =  0.  For 
a  smooth  loading  pulse  the  nominal  axial  strain  rate  and  strain  acceleration  are  initially  zero  also,  4(0)  =  0  and 
es(0)  =  0;  and  4(i )  increases  smoothly  with  t  up  to  time  4 ,  after  which  4  is  remains  constant  at  the  test  strain  rate 
4i  =  4(4)  >  0.5  At  early  times,  the  strain  acceleration  term  in  (25),  e2,  dominates  as  4  increases  from  its  initial 
value  of  zero;  ez  eventually  reaches  a  peak  (positive)  value  and  then  decays  to  zero  as  4  approaches  the  plateau 
strain  rate  4:1-  This  peak  in  ez  results  in  a  corresponding  early  peak  in  the  radial  stress.  However,  while  the  radial 
stress  initially  decreases  after  this  peak,  it  does  not  decay  to  zero  with  ez  since  the  (4)2  term  in  (25)  is  positive. 
In  fact,  since  ez  is  increasing,  the  (1  —  e2)3  term  in  the  denominator  is  decreasing,  and  hence  the  3(4)2/(l  —  e2)3 
term  is  strictly  increasing,  even  in  the  plateau  region  when  4(i )  =  4:1  ■  Consequently,  orr(t)  begins  to  increase  just 
prior  to  time  4  and  continues  to  do  so  while  4  remains  constant,  i.e. ,  there  is  a  strain  amplification  effect  on  the 
inertially  generated  radial  stress.  These  features  are  illustrated  in  Figure  1,  where  the  rise  time  4  =  135  ps  and 
the  plateau  value  of  the  nominal  axial  strain  rate  is  4i  =  2500/s.  For  this  strain  rate  history,  e2(4)  =  0.17  and 
e2(300^is)  =  0.58.  Note  that  by  (25)i,  the  product  Ar-A^  is  the  mean  radial  stress  arr  normalized  by  the  factor 
Po Rq  /4;  it  is  this  normalized  mean  radial  stress  that  is  plotted  in  the  figure.  Also  note  that  by  (25)i  and  (20)2,  the 
radial  stretch  acceleration  Ar  is  the  mean  value  of  the  nominal  radial  stress  Er#  normalized  by  the  same  factor. 


For  sufficiently  soft  materials  and  sufficiently  high  strain  rates,  the  inertial  effects  discussed  above  must  be  taken 
into  account  when  analyzing  the  data  from  SHPB  tests.  In  this  regard,  the  following  relations  for  the  mean  values 
of  the  stress  components  are  useful: 

O  zz  =  §  szz  +  &rr  i  Szz  =  ^  (tTzz  —  (7rr)  ,  <7 22  =  (1  —  e2)  E Zz  ■  (26) 

4However,  reducing  the  specimen  radius  also  reduces  the  signal  to  the  transmission  bar,  so  substantial  reduction  in  specimen  size 
must  be  accompanied  by  a  corresponding  reduction  in  the  diameter  of  the  pressure  bars.  This  is  one  of  the  motivations  for  the  use  of 
miniaturized  Hopkinson  bars. 

5For  simplicity,  we  neglect  the  subsequent  unloading  stage  in  both  the  analysis  and  the  numerical  simulations. 
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Figure  2:  Axial  strain  and  radial  stretch  rates  (a)  and  accelerations  (b)  for  a  smooth  ramp-up  to  a  constant  nominal 
radial  stretch  rate  of  1747/s  after  135  p, s.  Also  shown  in  (b)  is  the  normalized  mean  radial  stress  (dashed  line). 


These  follow  from  (14)  and  (20)i,  and  may  be  regarded  either  as  volumetric  averages  or  as  cross-sectional  averages; 
the  two  are  equivalent  if,  as  will  be  assumed  here,  the  specimen  is  in  dynamic  equilibrium  so  that  the  axial  stresses 
are  axially  uniform.  The  well-known  relation  on  the  right  expresses  the  mean  axial  Cauchy  stress  in  terms  of  the 
axial  strain  and  the  mean  axial  nominal  stress,  both  of  which  are  measured  (or  inferred  from  measurements)  in  an 
SHPB  test.  The  relation  on  the  left  in  (26)  implies  that  if  the  early  spike  in  arr  is  sufficiently  large  relative  to  szz 
(a  situation  that  could  occur  for  sufficiently  high  strain  rates  and  sufficiently  soft  specimens),  then  a  corresponding 
spike  in  the  measured  axial  stress  is  to  be  expected.  These  inertial  spikes  have  indeed  been  observed  in  SHPB  tests 
on  soft,  nearly  incompressible  materials;  cf.  [4],  [8],  [9],  and  [10].  The  middle  relation  in  (26),  which  is  equivalent  to 
the  relation  on  the  left,  indicates  that  an  “inertial  correction”  must  be  applied  to  the  (quasi-static)  uniaxial  stress 
relation  szz  =  |  azz  when  arr  is  sufficiently  large.  Since  arr  can  be  determined  from  the  measured  axial  strain  rate  4 
by  (25)3,  the  relation  (26)2  could  provide  a  means  to  estimate  the  axial  component  of  the  deviatoric  stress,6  provided 
the  assumptions  on  which  our  analysis  is  based  are  approximately  valid  and  the  difference  between  azz  and  dyr  is 
not  so  small  that  it  is  in  the  noise  level.  Then  the  radial  and  hoop  components  of  deviatoric  stress  can  be  determined 
from  (13)i. 

4.  Optimal  Loading  Pulses 

Since  a  constant  nominal  axial  strain  rate  does  not  eliminate  inertially  generated  radial  and  hoop  stresses,  it  is 
reasonable  to  seek  axial  strain  rate  histories  that  do  so.  On  setting  arr(t)  =  0  in  (25),  we  see  that  the  expression 
on  the  right  yields  an  ODE  for  ez:  2ez  +  3(4)2/(l  —  ez)  =  0.  However,  it  is  simpler  to  proceed  as  follows.  From 

(25)i  we  see  that  arr  =  0  iff  Xr  =  0  iff  Xr  is  constant  iff  er  is  constant.  Since  the  loading  pulse  arrives  at  t  =  0  and 
since  Ar  =  1  in  the  undeformed  state,  we  take  Ar(0)  =  1.  For  a  smooth  loading  pulse  we  cannot  impose  the  constant 
radial  strain  rate  condition  initially.  Instead  we  want  Ar(0)  =  0  and  Ar(0)  =  0,  with  Xr(t)  increasing  smoothly  and 
monotonically  with  t  up  to  some  time  t\ .  after  which  Ar  remains  constant: 

A r(t)  =  Xr(ti)  >  0 ,  for  t  >  ti .  (27) 

Then  Az(f)  =  1/A/  [cf.  (19)2],  and  from  A2  we  can  determine  ez,  4  and  e2;  alternatively,  we  can  determine  4  and 
ez  directly  from  Ar  and  its  rates  by  using  (24)  and  the  middle  relation  on  (21).  Since  A r(t)  is  positive  for  t  >  0,  so 
is  ez(t)  [cf.  (21)].  However,  A r(t)  =  0  for  t  >  t1;  so  by  (24)  we  see  that  ez(t)  <  0  for  t  >  hence  ez(t)  is  decreasing 
for  t  >  t\.  Since  the  initial  condition  Ar(0)  =  0  implies  4(0)  =  0,  it  follows  that  4(t)  must  increase  from  zero  to  a 
peak  value  at  some  time  tp  <  t\ ,  after  which  4  decreases.  In  spite  of  this  non-standard  feature,  the  resulting  nominal 
axial  strain  rate  history  ez{t)  will  have  the  property  that  for  t  >  ti,  arr(t)  =  0  and  hence  [cf.  (17)]  arr(r,z,t )  =  0 
throughout  the  specimen. 

®This  is  essentially  what  was  done  in  Sanborn  [10]  for  SHPB  tests  on  various  rubbers.  His  “corrected”  axial  stress  is  azz  —  arr  and 
hence  is  equivalent  to  the  |sz2  term  in  (26)  1 . 
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Figure  3:  Mean  values  of  the  axial,  radial  and  hoop  stress  and  of  3/2  the  axial  deviatoric  stress,  for  the  axial  strain 
rate  history  in  Figure  1.  The  radial  and  hoop  stresses  are  indistinguishable. 

For  comparison  with  the  more  conventional  strain  rate  history  considered  in  Figure  1,  we  take  the  rise  time  to 
the  constant  radial  strain  rate  to  be  the  same  as  the  previous  rise  time  to  the  constant  axial  strain  rate,  namely 
tl  =  135  /x s;  and  we  chose  the  plateau  value  Ar(ti)  =  1747/s,  since  this  results  in  ez(t\)  =  2500/s.  The  resulting  axial 
strain  and  radial  stretch  rates  are  shown  in  Figure  2. a.  The  corresponding  strain  and  stretch  acceleration  histories 
are  shown  in  Figure  2.b  along  with  the  normalized  mean  radial  stress,  which  is  nearly  indistinguishable  from  Ar.  For 
this  strain  rate  history,  ez(t\)  =  0.20  and  ez(350/ts)  =  0.55. 

9.  Numerical  Simulations 

To  test  the  optimal  pulse  shaping  and  the  inertial  correction  theory,  we  performed  numerical  simulations  of  hypo¬ 
thetical  SHPB  tests  on  a  soft,  nearly  incompressible,  solid  specimen  using  the  Lagrangean,  3-D  finite  element  code 
PRESTO  from  Sandia  Laboratories.  The  initial  radius  of  the  specimen  was  Rq  =  6.35  mm  and  the  initial  thickness 
was  Lq  =  1.45  mm,  giving  a  lengtli-to-diameter  ratio  of  0.11.  The  incident  and  transmission  bars  were  included 
in  the  simulation,  and  the  specimen-bar  interfaces  were  treated  as  frictionless.  We  used  1000  mm  long  aluminum 
bars  with  a  radius  of  19.05  mm.  The  specimen  and  bar  dimensions  (except  for  the  bar  lengths)  were  taken  from  the 
experimental  study  on  gelatins  by  Moy  et  al.  [8].  An  isotropic,  nonlinear  elastic  model  was  used  for  the  specimen 
(a  compressible  version  of  the  Mooney-Rivlin  model).  The  model  was  calibrated  to  give  rough  agreement  with  the 
large  strain,  quasi-static,  uniaxial  compression  data  on  the  20%  ballistic  gelatin  tested  in  [8]. 7  This  resulted  in  a 
small  strain  shear  modulus  of  80  kPa.  Since  20%  ballistic  gelatin  is  80%  water,  we  used  the  bulk  modulus  of  water, 
2.3  GPa;  the  ratio  of  bulk  to  shear  modulus  is  2.9  x  104. 

The  loading  wave  was  generated  by  an  imposed  axial  velocity  history  vz{t)  at  the  far  end  of  the  incident  bar.  This 
velocity  was  obtained  from  the  desired  nominal  axial  strain  rate  4,  using  the  approximate  relation 

vz{t)  «  ±L0ez(t) ,  (28) 

which  neglects  the  motion  of  the  specimen-transmission  bar  interface  and  assumes  the  particle  velocity  doubles  on 
reflection  from  the  specimen-incident  bar  interface.  The  nominal  axial  strain  rate  histories  used  in  computing  vz{t) 
were  those  from  Figures  1  and  2,  but  since  the  relation  (28)  is  only  approximate,  the  actual  strain  rates  and  strain 
accelerations  in  the  specimen  would  differ  somewhat  from  those  in  the  figures  even  if  the  specimen  deformation  were 
approximately  homogeneous. 

Figures  3  and  4  plot  the  histories  of  azz ,  arr,  agg  and  | szz.  These  are  the  (volumetric)  mean  values  of  the  radial, 

7Ballistic  gelatin  is  a  viscoelastic  material,  and  the  quasi-static  tests  showed  some  increase  in  stress  with  increasing  strain  rate.  The 
model  was  calibrated  so  that  the  axial  stress-strain  curve  was  slightly  above  the  curve  for  highest  (quasi-static)  strain  rate  of  1/s. 


Figure  4:  Mean  values  of  the  axial,  radial  and  hoop  stress  and  of  3/2  the  axial  deviatoric  stress,  for  the  axial  strain 
rate  history  in  Figure  2.  The  radial  and  hoop  stresses  are  indistinguishable. 


hoop  and  axial  stress  and  3/2  the  mean  value  of  the  axial  deviatoric  stress,  as  computed  in  the  simulations.  Recall 
that  azz  =  |  szz  for  a  uniaxial  stress  state;  cf.  (26) i  with  arr  set  to  zero.  The  times  in  Figures  3  and  4  are  shifted 
relative  to  those  in  Figures  1  and  2,  since  in  the  simulations  t  =  0  is  the  instant  at  which  the  velocity  is  applied  at 
the  far  end  of  the  incident  bar.  Figure  3  is  for  a  nominal  axial  strain  rate  history  given  (approximately)  by  that 
in  Figure  l.a,  that  is,  when  e,  is  eventually  constant.  The  radial  and  hoop  stresses  are  indistinguishable,  and  they 
increase  substantially  after  the  inertial  spike.  The  compressed  specimen  began  squeezing  out  beyond  the  bars  at 
around  453  /is. 

Figure  4  is  for  a  nominal  axial  strain  rate  history  given  (approximately)  by  that  in  Figure  2. a,  that  is,  when  Ar  is 
eventually  constant  (an  optimal  pulse  shape).  The  radial  and  hoop  stresses  are  again  indistinguishable,  but  now 
they  drop  to  nearly  zero  after  the  inertial  spike,  that  is,  once  Xr  reaches  its  plateau  value,  which  occurs  at  a  nominal 
strain  of  0.20.  Thereafter,  azz  is  very  close  to  |  szz.  Both  of  these  facts  indicate  that  a  nearly  uniaxial  stress  state 
has  been  achieved.  When  comparing  Figures  3  and  4,  keep  in  mind  that  the  scales  on  the  vertical  axis  are  different. 
The  inertial  spike  in  Figure  4  is  only  slightly  larger  than  that  in  Figure  3.  For  the  case  considered  in  Figures  1  and  3, 
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Figure  5:  The  computed  mean  axial  stress  (szz)  compared  with  two  estimates  for  it. 

Figure  5  provides  a  check  on  the  inertial  correction  (26)2  for  the  axial  deviatoric  stress:  szz  =  |  (crzz  —  arr).  In  the 


figure  legend,  szz  denotes  the  mean  axial  deviatoric  stress  szz  as  computed  in  the  simulation,  s22|  uniaxial  denotes  an 
estimate  for  the  deviatoric  stress  that  is  only  valid  for  a  uniaxial  stress  state,  namely  |  azz.  It  clearly  disagrees  with 
szz  and  shows  the  error  that  would  result  in  neglecting  inertial  effects  in  the  analysis  of  the  data.  Finally,  s22|corrected 
denotes  the  estimate  for  the  deviatoric  stress  obtained  from  the  inertial  correction  above,  using  azz  computed  in  the 
simulation  and  arr  determined  from  (25)3.  It  is  quite  close  to  the  actual  mean  value  for  most  of  the  simulation. 
Near  the  end  of  the  simulation  the  conditions  departed  substantially  from  those  on  which  the  analysis  was  based: 
the  specimen  bulged  and  was  no  longer  in  dynamic  equilibrium. 

10.  Discussion  and  Conclusions 

The  results  of  the  analysis  and  the  preliminary  numerical  simulations  indicate  that  radial  inertial  effects  in  SHPB 
tests  on  soft,  nearly  incompressible  materials  could  be  eliminated  (after  a  preliminary  inertial  spike)  by  tailoring 
the  loading  pulse  so  that  the  nominal  radial  strain  rate  in  the  specimen  is  eventually  constant.  The  corresponding 
axial  strain  rate  history  is  easily  determined  analytically  and  can  be  easily  (though  only  approximately)  imposed 
computationally.  Whether  such  loading  pulses  can  be  generated  experimentally  with  use  of  pulse  shapers  remains 
to  be  seen.  Recently,  Casern  [14]  has  developed  a  technique  for  tailoring  loading  pulse  shapes  by  means  of  graded 
impedance  striker  bars.  This  method,  in  conjunction  with  conventional  pulse  shapers,  may  possibly  provide  a  means 
to  generate  the  pulse  shapes  considered  here. 
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